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1.0 


INTRODUCTION 


One  is  familiar  with  the  display  of  a  passive  sonar  signal 
as  a  lofargram,  where  the  "instantaneous"  spectrum  is  displayed 
horizontally  as  a  function  of  time  running  vertically  downward.  Of 
course,  there  is  actually  no  such  thing  as  an  instantaneous  spec¬ 
trum;  rather,  the  spectrum  is  determined  for  a  time  interval  At  , 
either  by  taking  the  discrete  Fourier  transform  of  the  signal 
sampled  during  that  time  Interval  at  the  appropriate  Nyquist  rate, 
or  else  by  analog  means  using  narrowband  filters. 

Also  generally  known  is  that  the  frequency  resolution  of 
the  spectrum.  At  ,  is  not  independently  specifiable,  but  rather 
varies  inversely  to  the  time  interval  according  to  (what  physicists 
call)  the  uncertainty  principle. 

At  Af  -  1  .  (1.1) 

Thus,  one  can  choose  to  have  either  lofargrams  with  fine  time 
resolution  and  poor  frequency  resolution  (obtained  using  short 
coherent  integration  times)  or  else,  vice  versa,  with  fine  frequency 
resolution  and  poor  time  resolution  (obtained  using  long  coherent 
integration  times). 


The  purpose  of  this  note  is  to  provide  a  pedagogical 
introduction  to  the  so-called  "Wigner  distribution  function" ^  ^ ,  an 
to  demonstrate  that  what  is  generally  known  is  not  necessarily 
correct:  the  Wigner  distribution  function  displays  a  signal  as  a 
function  of  both  frequency  and  time,  but  in  a  manner  that  is  not 
resolution  limited  in  either  domain.  We  think  it  possible  that  the 
Wigner  distribution  can  be  used  to  provide  a  graphic  realization  of 
passive  sonar  signals  that  is  sufficiently  close  in  character  to 
conventional  lofargrams  so  that  their  body  of  accumulated  wisdom 
will  still  be  applicable,  yet  provide  significant  new  detailed 
information  on  the  sonar  signal. 

Originally  invented  by  Wigner  and  Szilard^  for  appli¬ 
cations  in  quantum  statistical  mechanics  (see,  e.g.,  [2]),  the 
Wigner  distribution  has  surfaced  from  time  to  time  in  fluid 
dynamical  applications  (e.g.,  [3]),  in  optical  signal  processing 
[4,5],  and  for  the  analysis  of  speech  [8].  As  we  will  see,  it  is 
closely  related  to,  but  different  from,  some  standard  concepts  in 
radar  processing  [6]. 


2.0 


FOURIER  TRANSFORM  OF  THE  SPECTROGRAM  PLANE 


He  can  define  Che  iofargram  of  a  signal  s(t)  as  follows: 
Centered  around  a  given  time  t  »  multiply  the  signal  by  a 
localized  window  function  w( t  -  tQ)  ,  then  Fourier  transform  the 
result.  We  obtain  a  function  of  both  frequency  f  and  time  t  . 


s(f,to> 


/  dt  e2nift  s(t)w(t  -tj 


(2. 


The  window  function  w  is  frequently  taken  to  be  a  square  window  of 
some  width  t  ,  equivalent  to  taking  a  uniform  sample  of  the 
signal.  Alternatively,  the  window  might  be  a  Gaussian  of  width 
t  .  We  will  assume  this  latter  case  as  a  mathematical 
convenience . 


Something  Interesting  will  emerge  if  we  take  the  two- 
dimensional  Fourier  transform  of  (2.1) 


2w  iX  t  2u  iy  f 

R(A,u)-/dtoe  7  df  e  S(f,to)  (2. 


2»Ut  2ffi„f  2* ift  -2ir ift' 
"  ////  dt  df  dt  dtf  e  e  e  e 


S(t)8(t*)w(t  -  tQ)w(t'  -  tQ)  . 


The  second  line  of  (2.2)  follows  from  writing  out  the  modulus  square 
of  (2.1)  as  the  product  of  the  integral  with  its  complex  conju¬ 
gate.  The  signal  s(t)  is  assumed  to  be  real. 


The  integrals  over  f  and  t'  in  (2.2)  can  be  done,  giving  in 


turn 


R(X,p) 


2*iXt 

///  dtQ  dt  dt'  6(y+t-t')e  Os(t)s(t,)w(t-to)w(t,-to) 

(2.3) 


2idXt 

/  dt  s(t)s(t+p)  /  dtQe  °  w(t-tQ)w(t+(rt0)  . 


Now  notice  that  the  second  integral  in  the  second  line  of  (2.3) 
depends  only  on  the  window  function,  not  on  the  signal.  For  a 
Gaussian  window 

w(x)  -  e“*  (2.4) 


the  Integral  can  be  done  analytically,  giving 


R(  X,  p)  -  e 


-y2/(**2)„-( 


J  dt  s(t  -  -£>(t  +  -|)e 


2wLXt 


(2.5) 
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In  (2.5)  we  have  also  redefined  the  dummy  variable  t  so  as  to  make 
the  Integrand  more  symmetrical. 

Notice  that  the  width  x  of  the  window  function  appears 
only  in  the  exponentials  in  front  of  the  integral,  and  that  the 
integral  is  a  functional  only  of  the  signal  s(t)  and  of  the  para¬ 
meters  X  and  ii  .  If  we  make  t  very  large,  then  the 
exponentials  in  front  cut  off  very  rapidly  in  A  ,  but  only  very 
slowly  in  ]t  .  If  we  make  t  very  small,  the  reverse  is  true.  In 
general,  the  only  function  of  the  exponentials  in  front  of  the 
integral  is  to  "enforce"  the  uncertainty  principle  in  the  2- 
dimensional  Fourier  domain:  the  lofargram  cannot  simultaneously 
contain  rapidly  varying  information  in  the  f  (conjugate  variable 
p  )  and  tQ  (conjugate  variable  X  )  directions. 

But  why  not  simply  strip  off  the  exponentials  that  enforce 
the  uncertainty  principle  and  examine  the  remaining  piece,  which 
depends  only  on  the  signal, 

Q(A,p)  =  /  dt  s (t  -  |)s(t  +  ■£)  e2,rl^t  .  (2.6) 


In  fact,  this  object  has  a  conventional  name ;  it  is  the  symmetrical 
autoambiguity  function  of  the  signal  (see  Ref.  6).  So  the  (purely 
pedagogical)  result  of  this  section  is:  the  autoambiguity  function 
of  a  signal  is  the  (window-independent  piece  of  the)  2-dimensional 
Fourier  transform  of  its  lofargram. 


3.0 


WIGNER  DISTRIBUTION  FUNCTION 


Now  we  do  Che  obvious  thing.  We  Cake  inverse  Fourier 
transforms  Co  get  back  Co  Che  2-dimenslonal  lofargram  space.  We 
will  not  get  back  our  original  lofargram,  because  we  have  stripped 
off  Che  dependence  on  Che  window  function.  We  want  Co  emphasize 
Chat  this  is  not  Che  same  as  Caking  any  particular  window  (e.g., 
very  wide  or  very  narrow) ;  1C  corresponds  instead  Co  taking  no 
window  at  all! 

-2itit  X  -2irifp 

P(f,tQ)  =  I  e  °  d X  /  e  du  Q(X,y) 

(3.1) 

-/  du.(to-f).(to+f).2rtf“  . 

This  is  Che  Wlgner  distribution  function  of  Che  signal  s(t).  It  is 
a  function  of  both  time  and  frequency,  just  like  the  lofargram.  The 
second  line  of  (3.1)  is  obtained  by  doing  the  integral  over  X  in 
the  first  line. 

An  interesting  identity  is  that  the  Wigner  distribution 
P(f,tQ)  can  also  be  expressed  in  terms  of  the  Fourier  transform 
S(t)  of  the  signal  s(t) 


*  1  1  “2lrlXtn 

p(f.t0)  -  /S(f-|)S(f+|)e  dx  .  (3.2) 

We  see  that  the  Wigner  distribution  function  is  quadratic 
in  the  signal  amplitude,  just  like  the  lofargram  (2.1).  It  is  not, 
however,  positive  definite.  We  will  see  why  this  is  so  in 
section  4. 

It  is  worth  noting  that  the  computation  of  the  Wigner 
distribution  is  not  significantly  slower  than  the  computation  of  the 
lofargram  (2.1),  which  requires  one  Fourier  transform  (FFT)  for  each 
time  t  of  interest.  The  Wigner  distribution  (3.1)  also  requires 
one  FFT  for  each  value  of  tQ  .  The  series  of  points  that  is  FFT'd 
is,  however,  not  the  windowed  signal,  but  is  rather  the  symmetrical 
product  of  the  signal  at  positive  and  negative  lags. 

References  [4]  and  [5]  prove  many  theorems  about  the  Wigner 
distribution,  and  give  further  references.  We  will  not  repeat  that 
material  here.  Rather,  in  the  next  section,  we  will  sketch  a  couple 
of  simple  analytic  examples  which  give  some  Insight  into  what  one 
might  actually  "see"  in  the  Wigner  distribution  of  a  passive  sonar 
signal. 


4.0 


EXAMPLES 


4.1 


i 


Segment  of  Pure  Tone 


Suppose  the  input  signal  is 


s(t)  ■  cos(ut) 

-  0 


-i  <‘<i 


otherwise 


(4.1 


Then  the  Higner  distribution  P(f ,tQ )  is  evidently  zero  outside  of 
the  interval  -  4  <  t  <  —  ,  since  in  the  integral  (3.1)  one  or  the 
other  of  the  two  terms  in  s(t  +  yj  will  be  zero.  So,  we  see  that 
the  Wlgner  distribution  turns  on  Instantaneously  with  the  signal. 


In  the  interval  -  4  <  t  <4  where  it  is  nonzero, 

2  0  2 

P(f,tQ)  is  symmetric  around  t  ■  0,  so  without  loss  of  generality 
we  can  confine  attention  to  the  case  of  0  <  t  <  4  .  In  this  case 

O  L 

U  ranges  from  -  (  y  -  t  )  to  +  (  y  -  t  )  .  Using 
trigonometric  identities,  (3.1)  can  be  rewritten  as 


dy  cos  (2vf)i)  [cos(2u»to)  +  cos(u>ii)] 


sin[(u+2wf ) (■j  -  tQ)]  sin[(u-2wf)(-j  -  tQ)] 


u  +  2ir£ 


-  2»f 


Notice  that  (4.2)  is  the  sum  of  three  resonance  integrals.  Each  has 
a  main-lobe  width  in  the  frequency  direction  of  about  1/t  .  In 
other  words,  the  longer  the  duration  of  the  tone,  the  narrower  the 
frequency  resolution  reflected  in  the  Wigner  distribution.  This  is 
just  as  we  would  expect  it  to  be.  with  the  novel  feature  that  there 
is  no  need  for  us  to  search  for  the  "best"  integration  time  for  the 
signal:  the  Wigner  distribution  does  this  for  us  automatically. 


Where,  in  the  frequency  direction,  are  the  three 
resonances?  Two  are  at  f  -  +  u/2v  ,  where  we  would  expect  them  to 
be.  The  third  is  at  f  ■  0.  To  understand  this  one,  we  need  to 
examine  the  dependence  of  the  three  resonances  on  time  tQ  .  The 
dependence  on  tQ  inside  the  resonance  terms  just  gives  all  three 
resonances  a  kind  of  broadening  as  the  boundary  in  time  is 


s 


approached.  More  Important  is  the  rapidly  oscillating  term,  with 
angular  frequency  twice  that  of  the  signal,  which  multiplies  the 
zero-frequency  resonance  term.  Figure  1  sketches  the  location  and 
sign  of  the  main-lobe  peaks  or  ridges  of  equation  (4.2). 

To  understand  what  is  going  on  in  the  figure,  imagine  that 
the  figure  is  averaged  or  smeared  in  the  horizontal  direction.  Then 
the  amplitudes  at  zero  frequency  average  to  zero,  and  one  will  be 
left  with  narrow  frequency  resonances  at  f  »  +  <a/2n  which  are 
smeared  out  in  time.  This  is  precisely  the  lofargram  of  the  signal 
if  constructed  with  a  wide  time-window:  it  maintains  frequency 
resolution  but  smears  out  time  resolution.  By  contrast,  next 
imagine  that  the  figure  is  averaged  or  smeared  in  the  vertical 
direction.  Then  the  zero-frequency  positive  peaks  reinforce  the 
positive  values  at  f  -  +  vjlit  to  give  vertical  smears  located  at 
the  positive  and  negative  peaks  of  the  original  signal,  while  the 
zero-frequency  negative  peaks  exactly  cancel  with  the  positive 
values  at  f  •  +  <V2*  to  give  zero  response  at  the  zeros  of  the 
original  signal.  This  is  exactly  the  lofargram  of  the  signal  if 


constructed  with  a  narrow  time-window 


So,  in  summary,  we  see  that  the  Wigner  distribution 
displays  simultaneously  all  the  structures  which  would  emerge  from 


lofargram  with  any  window  width,  or  (for 
window  shape. 


that  matter)  with  any 
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.  <V'. 


4.2 


Periodic  Series  of  Clicks 


We  can  see  the  saae  general  ideas  emerging  in  another 
example.  Consider  the  signal 


+00 

s(t)  -  l  6(t  -  nT)  . 
n*-n 


(4.3) 


This  represents  a  sequence  of  sharp  clicks  (delta  functions), 
equally  spaced  with  period  T.  The  Wlgner  distribution  is  thus 


JdU  l  «(t0  -«»T)  «(t0+|-nT)e2l,lfiJ 

n 

-  I  e2*if<2nT>  for  t0  -  ...-2T,-T,0,T,2T... 

n 

y  e2xif(2n+l)T  for  t  -  2  T,- -  T,~  T,- T... 

n  °  2  2  2  2 

■  l  I  «(v  «)  «(*  -  2f)  +  <_l>  ‘  ‘(v  lk  +  li1) 
k  ‘  (»..) 


(except  for  a  normalizing  constant).  Pigure  2  shows  the  locations 
and  signs  of  the  delta-function  peaks  out  of  which  this  Wlgner 
distribution  is  composed. 


Notice  that  the  positive  delta  functions  form  two  inter¬ 
locking  sets  of  lines,  horizontal  and  vertical.  In  the  centers  of 
the  lattice  thus  formed  are  the  Isolated  negative  delta  functions. 
The  vertical  lines  of  positive  delta  functions  are  separated  in  time 
by  T,  and  they  represent  the  signal  with  infinitely  good  time 
resolution.  The  horizontal  lines  of  delta  functions  are  separated 
in  frequency  by  a  spacing  of  1/T,  and  they  represent  the  fundamental 
and  higher  harmonics  of  the  signal  (which  is,  after  all,  perfectly 
periodic)  with  infinitely  good  frequency  resolution. 

So  what  function  do  the  negative  delta  functions  perform? 
The  answer  can  be  seen  by  again  imagining  that  the  Wigner  function 
is  smeared  either  horizontally  or  vertically.  In  that  case,  the 
negative  values  cancel  out  either  the  power  between  the  perfect 
harmonic  rows,  or  the  power  between  the  perfectly  time-resolved 
pulses,  in  both  cases  giving  realizable  lofargrams. 

In  fact,  one  can  see  intuitively  that  any  smearing  of 
Figure  2  that  obeys  the  uncertainty  principle  will  give  a  positive- 
definite  lofargram:  Since  each  minus  sign  is  surrounded  by  eight 
plus  signs,  there  is  no  way  to  center  a  window  of  unit  area  on  the 
minus  sign  without  bringing  in  more  than  enough  plus  signs  to  give 
positivity.  It  is  generally  true  on  Wigner  distributions  that 


5.0 


DISCUSSION 


It  may  require  some  experimentation  to  discover  the  best 
graphical  display  of  a  Wigner  distribution  function.  We  have  not 
yet  done  this  experimentation,  but  intend  to  do  so  in  the  near 
future.  One  desires  a  format  in  which  the  information  is  presented 
in  a  manner  as  close  as  possible  to  that  of  the  standard  lofargram. 

One  good  candidate  for  a  display  format  is  a  grey-scale 
presentation  ranging  from  the  most  negative  values  (pure  white)  to 
most  positive  (pure  black).  In  this  format,  regions  where  negative 
and  positive  lobes  are  in  close  proximity  will  have  a  mid-scale  grey 
tone,  while  positive  features  will  appear  as  darker  lines  or 
regions.  The  human  eye  is  quite  good  at  picking  out  linear  black 
features  in  a  noisy  grey-scale  background,  so  that  details  of  high 
resolution  simultaneously  in  frequency  and  time  may  be  readily 
apparent . 


Other  display  formats  are  also  under  consideration  and  will 
be  tried  in  due  course. 

Although  we  do  not  yet  have  actual  data,  we  have  reason  to 
hope  for  a  graphical  presentation  in  which  unrelated  signals  will 


appear  tilth  distinguishable  natural  line  widths,  so  the  discrimi¬ 
nation  between  submarine  and  surface-ship  signatures  is  aided. 

Also,  we  would  hope  that  the  detectability  of  very  weak  but  very 
narrow  lines  will  be  furthered,  if  they  appear  "automatically"  in  a 
display,  without  the  necessity  of  special  narrowband  processing. 
Finally,  when  high-resolution  frequency  and  high-resolution  time 
information  is  displayed  simultaneously,  there  may  be  serendipitous 
features  in  the  signal  which  would  not  otherwise  be  noted,  and  which 
may  yield  new  discrimination  techniques. 
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